Photon-Hadron Jet Correlations in p-\-p and Au+Au Collisions at y/sNN=200 GeV 



A. Adare," S. Afanasiev,^^ C. Aidala,!^'^^ N.N. Ajitanand,^'^ Y. Akiba,48.49 h. Al-Bataineh,'*3 J. Alexander 
A. Al-Jamel,^3 k. Aoki.^i'^s Aphecetche,^^ R. Armendariz,43 S.H. Aronson,^ J. Asai,''^'''^ E.T. Atomssa, 
R. Averbeck,55 T.C. Awes,'^'' B. Azmoun,*^ V. Babintsev,^! M. Bai,'^ G. Baksay,!"^ L. Baksay/"^ A. Baldisseri,- 
K.N. Barish/ P.D. Bames,^* B. Bassalleck,'*^ A.T. Basye,^ S. Bathe/ S. Batsouli,^^''*'' V. Baublis,^^ F. Bauer/ 
C. Baumann/** A. Bazilevsky,^ S. Belikov/'^^^Q R. Bennett/^ A. Berdnikov/i Y. Berdnikov/i A.A. Bickley," 
M.T. Bjorndal/2 J.G. Boissevain/^ H. Borel/^ K. Boyle/^ M.L. Brooks/'^ D.S. Brown/^ D. Bucher/^ 
H. Buesching,6 V. Bumazhnov/i G. Buiice,*^'''^ J.M. Burward-Hoy/^ S. Butsyk/^.ss c.M. Camacho/" 
S. Campbell/5 J.-S. Chai/*^ B.S. Chang/^ w.C. Chang/ J.-L. Charvet/"* S. Chernichenko/i J. Chiba/^ 

C.Y. Chi/2 M. Chiu/2'22 I. J. Choi/3 R.K. Choudhury/ T. Chujo/^-eo P. Chung/^ A. Churyn/^ 
V. Cianciolo/'* Z. Citron/^ C.R. Cleven/^ Y. Cobigo/^ B.A. Cole/^ M.P. Comets/^ P. Constantin/^.s" 
M. Csanad/6 T. Csorgo/^ T. Dahms/^ S. Dairaku/i'^s k. Das/^ G. David/ M.B. Deaton/ K. Dehmelt/^ 

H. Delagrange/'^ A. Denisov/^ D. d'Enterria/^-^^ A. Deshpande/^' E.J. Desmond/ O. Dietzsch/^ A. Dion/^ 

M. DonadeUi/^ J.L. Drachenberg/ O. Drapier/^ A. Drees/^ K.A. Drees/ A.K. Dubey/^ A. Durum/^ 
D. Dutta/ V. Dzhordzhadze/^ " Y.V. Efremenko/'* J. Egdemir/^ F. EUinghaus/^ W.S. Emam/ T. Engelmore/^ 
A. Enokizono/o>33 H. En'yo/*'*^ B. Espagnon/^ S. Esumi/^ K.O. Eyser/ B. Fadem/^ D.E. Fields/^-^s 
M. Finger, Jr./' 25 M. Finger/' ^5 F. Fleuret/^ S.L. Fokin/" B. Forestier/^ Z. Fraenkel/^'Q J.E. Frantz/^'^s 
A. Franz/ A.D. Frawley/^ K. Fuji-wara/^ Y. Fukao/i''** S.-Y. Fung/ T. Fusayasu/^ S. Gadrat/^ 

I. Garishvih/^ F. Gastineau/^ M. Germain/^ A. Glenn/i'^^ H. Gong/^ M. Gonin/^ J. Gosset/" Y. Goto/*''*^ 
R. Granier de Cassagnac/^ N. Grau/2'24 S.V. Greene/° M. Crosse Perdekamp/^'^^ T. Cunji/" H.-A. Custafsson/^ 

T. Hachiya/O'-^s A. Hadj Henni/^ C. Haegemann/^ J.S. Haggerty/ M.N. Hagiwara/ H. Hamagaki/o R. Han/^ 
H. Harada/" E.P. Hartouni/^ K. Haruna/° M. Harvey/ E. Haslum/^ K. Hasuko/^ R. Hayano/° M. Heffner/^ 
T.K. Hemmick/5 T. Hester/ J.M. Heuser/^ X. He/^ H. Hiejima/^ J.C. Hill/* R. Hobbs/^ M. Hohlmann/^ 
M. Holmes/'' W. Holzmann,^* K. Homma/° B. Hong/^ T. Horaguchi/O'^s. 58 D. Hornback/^ S. Huang/" 
M.G. Hur/6 T. Ichihara/^' ^9 R. Ichimiya/^ Y. Ikeda/^ K. Imai/^'''* J. Imrek/5 M. Inaba/^ Y. Inoue/°'48 
D. Isenhower/ L. Isenhower/ M. Ishihara/^ T. Isobe/'^ M. Issah/'' A. Isupov/5 D. Ivanischev/'' B.V. Jacak/5'0 
J. Jia/2 J. Jin/2 o. Jinnouchi/9 B.M. Johnson/ K.S. Joo/° D. Jouan/5 F. Kajihara/°'''8 S. Kametani/O'^s.ei 

M TroTnilnOT-o 48,49, 58 T T^OTni,, 55 IV/r T^ono+o 49 T TT T^^CTin- 63 TT XToTirM, 48, 58 T Xt^orimoi 



Jia/^ J. Jin/'' O. Jinnouchi/9 B.M. Johnson/ K.S. Joo,"" D. Jouan/^ F. Kajihara/'^'^" S. Kametani/U'^'"^ 
Kamihara/8'49'58 J. Kamin/5 M. Kaneta/^ J.H. Kang/^ H. Kanou/8'58 J. Kapustinsky/" T. Kawagishi/^ 
D. Kawall/^'''^ A.V. Kazantsev/° S. Kelly/^ T. Kempel/'^ A. Khanzadeev/^ K.M. Kijima/" J. Kikuchi/i 
I. Kim/9 D H Kim/o D.J. Kim/^ E. Kim/^ S.H. Kim/^ Y.-S. Kim/^ E. Kinney/^ K. Kiriluk/i A. Kiss/^ 
Kistenev/ A. Kiyomichi/^ J. Klay/^ C. Klein-Boesing/^ L. Kochenda/"^ V. Kochetkov/^ B. Komkov/"^ 
, Konno/^ J. Koster/^ D. Kotchetkov/ A. Kozlov/^ A. Kral/^ A. Kravitz/^ P.J. Kroon/ J. Kubart/'^^ 
I J. Kunde/4 N. Kurihara/" K. Kurita/°'*8 ]y[ Kurosawa/*^ M.J. Kweon/^ Y. Kwon/'^'SS G.S. Kyle/^ 
3. Lai/2 Y.S. Lai/^ J.C. Lajoie/* D. Layton/^ A. Lebedev/'' Y. Le Bornec/5 S. Leckey/5 
Lee/9 M.K. Lee/^ T. Lee/^ M.J. Leitch/* M.A.L. Leite/^ B. Lenzi/^ P. Liebing/^ H. Lim,' 



D 
J 

N 

D. 

B. I. Kim/9 D H Kim/o q j Kim/^ E. Kim/^ S.H. Kim/^ Y.-S. Kim/^ E. Kinney/^ K. Kiriluk/i A. Kiss/^ 
E. Kistenev/ A. Kiyomichi/^ J. Klay/^ C. Klein-Boesing/^ L. Kochenda/"^ V. Kochetkov/^ B. Komkov/"^ 
M. Konno/9 J. K 

G.J. Kunde/4 N 

D.M. Lee/* K.B. Lee/^ M.K. Lee/^ T. Lee/^ M.J. Leitch/* M.A.L. Leite/^ B. Lenzi/^ P. Liebing/9 H. Lim/^ 
T. Liska/3 A. Litvinenko/5 H. Liu/^ M.X. Liu/* X. Li/ X.H. Li/ B. Love/" D. Lynch/ C.F. Maguire/° 
Y.L Makdisi/'*^ A. Malakhov/5 M.D. Malik/^ V.L Manko/" E. Manuel/^ Y. Mao/^'*^ L. Masek/'^s H. Masui/^ 
F. Matathias/2' 55 m.C. McCain/^ M. McCumber/5 P.L. McCaughey/* N. Means/5 B. Meredith/^ Y. Miake/^ 
P. Mikes/' 23 K. Miki/9 T.E. Miller/" A. Milov/'55 S. Mioduszewski/ G.C. Mishra/9 M. Mishra/ J.T. Mitchell/ 
M. Mitrovski/* A.K. Mohanty/ Y. Morino/° A. Morreale/ D.P. Morrison/ J.M. Moss/* T.V. Moukhanova/° 
D. Mukhopadhyay/o J. Murata/°'*^ S. Nagamiya/^ Y. Nagata/" J.L. Nagle/^ M. Naghs/2 M.L Nagy/^ 
L Nakagawa/8'*9 Y. Nakamiya/" T. Nakamura/" K. Nakano/8'58 J. Newby/^ M. Nguyen/5 T. Niita/'^ 
B.E. Norman/* R. Nouicer/ A.S. Nyanin/o J. Nystrand/^ E. O'Brien/ S.X. Oda/" C.A. Ogilvie/* H. Ohnishi/^ 
LD. Ojha/° H. Okada/i'*^ K. Okada/'' M. Oka/" 0.0. Omiwade/ Y. Onuki/^ A. Oskarsson/^ L Otterlund/*^ 
M. Ouchida/° K. Ozawa/° R. Pak/ D. Pal/" A.P.T. Palounek/* V. Pantuev/5 V. Papavassihou/^ J. Park/^ 
W.J. Park/9 g p Pate/3 H. Pei/* J.-C. Peng/2 jj Pereira/* V. Peresedov/5 D.Yu. Peressounko/" 

C. Pinkenburg/ R.P. Pisani/ M.L. Purschke/ A.K. Purwar/*'55 H. Qu/" J. Rak/*'*2 A. Rakotozafindrabe/2 

I. Ravinovich/2 K.F. Read/*'57 S. Rembeczki M. Reuter/5 K. Reygers/^ V. Riabov/^ Y. Riabov/'^ D. Roach/*^ 
C. Roche/5 s Rolnick/ A. Romana/2'B M. Rosati/* S.S.E. Rosendahl/^ P. Rosnet/5 P. Rukoyatkin/5 
P. Ruzicka/3 V.L. Rykov/^ S.S. Ryu/^ B. Sahlmueller/® N. Saito/^'**'*" T. Sakaguchi/' i°' S. Sakai/" 

K. Sakashita/8'58 H. Sakata/° V. Samsonov/^ H.D. Sato/i'*^ S. Sato/' 2^' 59 x. Sato/" S. Sawada/^ 
K. Sedgwick/ J. Seele/^ R. Seidl/2 A.Yu. Semenov/* V. Semenov/i R. Seto/ D. Sharma/2 T.K. Shea/ 
I. Shein/i A. Shevel/^'5* T.-A. Shibata/8'58 K. Shigaki/" M. Shimomura/9 T. Shohjoh/9 K. Shoji/i'*^ 



2 



P. Shukla,'' A. Sickles,^^'^^ c.L. Silva,^^ D. Silvermyr,^* C. Silvestre," K.S. Sim,^^ B.K. Singh,^ CP. Smgh,^ 
V. Singh,3 S. Skutnik,24 M. Slunecka,^^ ^5 q Smith,i A. Soldatov,^! R.A. Soltz,33 W.E. Sondheim,^^ 
S.P. Sorensen,57 I.V. Sourikova,'^ F. Staley,!'* P.W. Stankus ,"4 E. Stenlund,^^ M. Stepanov,''^ A. Ster,^^ 
S.P. Stoll,^ T. Sugitate,20 C. Suire,''^ A. Sukhanov,^ J.P. Sullivan,^^ J. Sziklai,^^ T. Tabaru,^^ S. Takagi,^^ 
E.M. Takagui,52 a. Taketani,48.49 r__ Tanabe,^^ j^.H. Tanaka,^^ Y. Tanaka,^! K. Tanida 48.49 m.J. Tannenbaum,*^ 
A. Taranenko,54 p. Tarjan/^ H. Themann,^^ T.L. Thomas,42 M. Togawa,3i.48 A. Toia,^^ J. Tojo 48 L. Tomasek,^^ 
Y. Tomita,59 H. Torii,2"'48 R.S. Towell/ V-N. Tram,32 1. Tserruya,62 Y. Tsuchimoto,20'48 S.K. Tuli,^ H. Tydesjo,^^ 
N. Tyurin,2i C. Vale,24 H. Valle,60 H.W. van Hecke,34 A. Veicht,^^ J. Velkovska,"^ R. Vertesi,!^ A.A. Vinogradov,3° 
M. Villus," V. Vrba,23 E. Vznuzdaev,47 M. Wagner,34'48 D. Walker,^^ X.R. Wang 43 Y. Watanabe,48.49 
F. Wei,24 J. Wessels,38 S.N. White,^ N. Willis,45 D. Winter,!^ C.L. Woody,^ M. Wysocki," W. Xie7'49 
Y.L. Yamaguchi,6i K. Yamaura,^^ R. Yang,22 A. Yanovich,^! Z. Yasin/ J. Ying,!^ S. Yokkaichi,48.49 q p.^ Young 44 
I. Younus42 I.E. Yushmanov,30 W.A. Zajc,^^ O. Zaudtke,^^ C. Zhang,i2.44 Zhou,^ J. Zimanyi,^^.^ and L. Zolin^s 

(PHENIX Collaboration) 

''Abilene Christian University, Abilene, TX 79699, U.S. 
^Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
'"^ Department of Physics, Banaras Hindu University, Varanasi 221005, India 
^'Bhahha Atomic Research Centre, Bombay 400 085, India 
^ Collider- Accelerator Department, Brookhaven National Laboratory, Upton, NY 11973-5000, U.S. 
''Physics Department, Brookhaven National Laboratory, Upton, NY 1 1973-5000, U.S. 
''University of California - Riverside, Riverside, CA 92521, U.S. 
^Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
^ China Institute of Atomic Energy (CIAE), Beijing, People's Republic of China 
Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

"University of Colorado, Boulder, CO 80309, U.S. 
^''Columbia University, New York, NY 10027 and Nevis Laboratories, Irvington, NY 10533, U.S. 
^'^ Czech Technical University, Zikova 4, 166 36 Prague 6, Czech Republic 
^'^Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
^^ELTE, Eotvos Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
"Florida Institute of Technology, Melbourne, FL 32901, U.S. 
Florida State University, Tallahassee, FL 32306, U.S. 
Georgia State University, Atlanta, CA 30303, U.S. 
^"Hiroshima University, Kagamiyama, Higashi- Hiroshima 739-8526, Japan 
IHEP Protvmo, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, 142281, Russia 
^'^ University of Illinois at Urbana- Champaign, Urbana, IL 61801, U.S. 
'' 'institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

^"''lowa State University, Ames, lA 50011, U.S. 
^'^ Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
'^^ KAERI, Cyclotron Application Laboratory, Seoul, Korea 
''^'KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
KFKI Research Institute for Particle and Nuclear Physics of the Hungarian Academy 
of Sciences (MTA KFKI RMKI), H-1525 Budapest 114, POBox 49, Budapest, Hungary 
^^Korea University, Seoul, 136-701, Korea 
'"''^Russian Research Center "Kurchatov Institute", Moscow, Russia 



31 



Kyoto University, Kyoto 606-8502, Japan 



"Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
'^'^ Lawrence Livermore National Laboratory, Livermore, CA 94550, U.S. 
'^'^Los Alamos National Laboratory, Los Alamos, NM 87545, U.S. 
^^LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
"^^ Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
^''Department of Physics, University of Massachusetts, Amherst, MA 01003-9337, U.S. 
^^Institut fiir Kernphysik, University of Muenster, D-48149 Muenster, Germany 
Muhlenberg College, Allentown, PA 18104-5586, U.S. 
Myongji University, Yongin, Kyonggido 449-728, Korea 
'^'Nagasaki Institute of Applied Science, Nagasaki-ski, Nagasaki 851-0193, Japan 
University of New Mexico, Albuquerque, NM 87131, U.S. 
^^New Mexico State University, Las Graces, NM 88003, U.S. 
^^Oak Ridge National Laboratory, Oak Ridge, TN 37831, U.S. 
^■'''IPN-Orsay, Universite Pans Sud, CNRS-IN2P3, BPl, F-91406, Orsay, France 
Peking University, Beijing, People's Republic of China 
^''PNPI, Petersburg Nuclear Physics Institute, Catchina, Leningrad region, 188300, Russia 



3 



^^RIKEN Nishina Center for Accelerator-Based Science, Wako, Saitama 351-0198, JAPAN 
^*"RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, U.S. 
^"Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 

'''''Saint Petersburg State Polytechnic University, St. Petersburg, Russia 
Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 
''''^System Electronics Laboratory, Seoul National University, Seoul, Korea 
Chemistry Department, Stony Brook University, Stony Brook, SUNY, NY 11794-3400, U.S. 
Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY 11794, U.S. 
SUBATECH (Ecole des Mines de Nantes, CNRS-IN2P3, Universite de Nantes) BP 20722 - 44307, Nantes, France 

^''University of Tennessee, Knoxville, TN 37996, U.S. 
Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
Vanderbilt University, Nashville, TN 37235, U.S. 
Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Ktkm-cho, Shmjuku-ku, Tokyo 162-0044, Japan 
^'^Wetzmann Institute, Rehovot 76100, Israel 
''^Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: March 19, 2009) 

We report the observation at the Relativistic Heavy Ion CoUider (RHIC) of suppression of back-to- 
back correlations in the direct photon-f-jet channel in Au-|-Au relative to p-\~p collisions. Two-particle 
correlations of direct photon triggers with associated hadrons are obtained by statistical subtraction 
of the decay photon-hadron (7-/1) background. The initial momentum of the away-side parton is 
tightly constrained, because the parton-photon pair exactly balance in momentum at leading order 
in perturbative quantum chromodynamics (pQCD), making such correlations a powerful probe of 
the in-medium parton energy loss. The away-side nuclear suppression factor, 7aa, in central Au-|-Au 
collisions, is 0.32 ± 0.12°*''* ± 0.09''^''* for hadrons of 3 < Px < 5 in coincidence with photons of 5 
< < 15 GeV/c. The suppression is comparable to that observed for high-px single hadrons 
and dihadrons. The direct photon associated yields in p-|-p collisions scale approximately with the 
momentum balance, zt = Pr / Pt ' expected for a measure of the away-side parton fragmentation 
function. We compare to Au-fAu collisions for which the momentum balance dependence of the 
nuclear modification should be sensitive to the path-length dependence of parton energy loss. 

PACS numbers: 13.85.Qk, 13.20.Fc, 13.20.He, 25.75.Dw 



I. INTRODUCTION 

Experimental results from RHIC have established the 
formation of hot and dense matter of a fundamentally 
new nature in relativistic heavy-ion collisions at ^/snn — 
200 GeV 1|. Energy loss in this dense nuclear mat- 
ter by color-charged, hard (E > 2 GeV) partons, and 
the jets into which they fragment, is generally accepted 
to be the mechanism responsible for the suppression of 
the high-pT hadron yields observed in central A-|-A col- 
lisions 0, [3- III the large multiplicity environment of 
heavy-ion collisions, two-particle correlations are often 
used to study jet modification and to infer properties of 
the medium. For example, high-px azimuthal dihadron 
correlations demonstrate that the degree of dijet away- 
side suppression depends on the pT of the "trigger" and 
"associated" hadrons. At moderate 3 GeV/c), the 

jet properties measured through two-particle correlations 
demonstrate novel features such as shape modifications 
which are thought to be a manifestation of the response 
of medium to the energy deposited by the attenuated 
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parton [1]. 

Di-hadron measurements of dijet pairs provide an am- 
biguous measurement of the energy loss of the away-side 
parton. The trigger hadron is a product of parton frag- 
mentation and therefore it is not possible to determine, 
event-by-event, whether the near-side parton has itself 
lost energy. Given the steeply falling jet spectrum, the 
sample of hard scatterings is biased towards configura- 
tions in which the parton loses little energy. In particu- 
lar, it is believed that hadron measurements are subject 
to a "surface bias" in which the hard scatterings sampled 
are likely to occur at the periphery of the overlap zone 
0, [B| ■ The away-side parton then is more likely to tra- 
verse a maximal path-length through the medium. For a 
sufficiently opaque medium, the attenuation of the par- 
ton may be nearly total, in which case the sensitivity to 
the average path-length is reduced Back-to-back, 
high-pT hadron pairs may originate preferentially from 
configurations in which the outgoing parton trajectories 
are tangential to the surface of the overlap zone Q. On 
the other hand, dihadron pairs may also originate from 
vertices deep in the collision zone if a parton has a fi- 
nite probability to "punch-through" or pass through the 
medium without interaction [9| . Calculations of the rela- 
tive importance of these two mechanisms depend both on 
the model of parton energy loss employed and the density 
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profile of the medium [1, [T^, [HI ■ 

Direct photon-jet pairs offer two major advantages in 
studying energy loss as compared to dijets because of 
the nature of the photon. First, in contrast to partons, 
photons do not carry color charge and hence do not in- 
teract strongly when traversing the medium p^ . The 
distribution of hard scattering vertices sampled by direct 
photon-triggered correlations is thus unbiased by the trig- 
ger condition. Suppression of the opposite jet is averaged 
over all path-lengths given by the distribution of hard 
scattering vertices. Second, at the Born level, direct pho- 
ton production in p+p and A-|-A collisions is dominated 
by the QCD Compton scattering process, q+g 9+7, 
and the photon momentum in the center-of-mass frame 
is exactly balanced by that of the recoil quark. Although 
higher order effects and other complications to this ideal- 
ized picture such as Next-to-Leading Order (NLO) 2 — > 3 
"fragmentation" photons or soft gluon radiation must be 
considered, the level of suppression can then be related 
directly to the energy loss of a parton of known initial 
momentum. In this way, the average path-length of the 
away-side parton may then be varied in a well controlled 
manner by selecting events of various momentum differ- 
ences between the j-h pair. 

For this reason, the 7+jet channel has long been con- 
sidered the "golden channel" for studying parton en- 
ergy loss [l^, [3] ■ Neglecting the above mentioned com- 
plications, specifically effects like transverse momentum 
broadening (the effect) and parton-to-photon frag- 
mentation, back-to-back j-h correlations in elementary 
collisions directly measure the fragmentation function of 
the recoil jet since z = p'^/pj"* « p^/p'^. In the standard 
picture of energy loss, partons are likely to lose some 
fraction of their energy in the medium, but are likely to 
fragment outside the medium. Hence, the parton energy 
loss can be considered an effective modification to the 
fragmentation function. Such a picture may be tested 
using "f-h correlations in nuclear collisions. Complemen- 
tary baseline measurements in p+p collisions are used to 
test the theoretical description of correlations in vacuum 
and to constrain possible contributions from higher or- 
der processes. Comprehensive reviews of direct photon 
phenomenology and data from elementary collisions may 
be found in [15,, J^, ,17,] . 



II. DETECTOR DESCRIPTION AND 
PARTICLE IDENTIFICATION 

The data were taken with the PHENIX detector [Hi 
using approximately 950 million Au+Au minimum bias 
events from the 2004 data set and 471 million photon- 
triggered events from the 2005 and 2006 p+p data sets 
corresponding to integrated luminosities of 3 (2005) and 
10.7 (2006) ph~\ The Beam-Beam Counters (BBC) [H, 
which are used to trigger the minimum bias data, select 
92% of the total inelastic cross section. In Au-|-Au the 
BBC and Zero-Degree Calorimeters (ZDC) were used for 



offline minimum bias event selection and centrality de- 
termination. In p+p collisions a high energy photon trig- 
ger, defined by coincidence between the BBC and a high 
energy Electromagnetic Calorimeter (EMCal) tower hit, 
was utilized. This EMCal based trigger [2(|| had an ef- 
ficiency of > 90% for events with photons and tt" with 
energies in the range used in the analysis and within the 
detector's geometric acceptance. 

The PHENIX central arms, each covering ±0.35 units 
of pseudorapidity around midrapidity and 90° in az- 
imuth, contain charged-particle tracking chambers and 
electromagnetic calorimeters [2l|. The EMCal [l^ con- 
sists of two types of detectors, six sectors of lead- 
scintillator (PbSc) sampling calorimeters and two of lead- 
glass (PbGl) Cerenkov calorimeters measuring EM en- 
ergy with intrinsic resolution as/E = 8.1%/y(i?)e2.1% 
and 5.9%/y/iE)® 0.8% respectively. The fine segmenta- 
tion of the EMCal (Ar? x A(/) ~ 0.01 x 0.01 for PbSc and 
- 0.008 X 0.008 for PbGl) allows for the reconstruction 
of 7r°'s and t^'s in the 27 decay channel out to pt of 20 
GeV/c. The details of direct photon, tt'~' and 77 meson 
detection and reconstruction within PHENIX have been 
described previously [I^ilHjH^I- Photon candidates with 
very high purity (> 98% for energies > 5 GeV) are se- 
lected from EMCal clusters with the use of cluster shower 
shape and charged particle veto cuts. Two-photon tt'^ and 
r] candidates are selected from photon pairs with pair in- 
variant mass in the appropriate tt" or ij mass range. Com- 
binatorial 27 background is reduced with cuts on energy 
asymmetry ai2 — \Ei — E2\/{Ei + E2), described in de- 
tail below. Some fraction of 7r° with starting at « 
13 GeV/c (in the PbSc detector) will appear as a single 
merged cluster, but with anomalous shower shape, and 
thus are removed from the analysis. The Tr'^and 77 mesons 
in the pT range from about 4 to 17 GeV/c and photons 
between 5 and 15 GeV/c are used in this analysis. For 
7 Pt between 13-15 GeV/c there is a < 2% contribu- 
tion of merged tt" cluster contamination, however this 
together with all sources of non-photon contamination 
are found to have a negligible impact on the two-particle 
correlation analysis of this report. Direct photons and 
their two-particle correlations are obtained by statistical 
subtraction of the estimated meson (mainly tt") decay 
photon contribution from the inclusive photon and 7-/1 
samples. 

Charged hadrons are detected with the PHENIX track- 
ing system 25| which employs a drift chamber in each 
arm spanning a radial distance of 2.0-2.4 m from the 
beam axis with a set of pixel pad chambers (PCI) di- 
rectly behind them. The momentum resolution was de- 
termined to be 5p/p — O.7%0 1.0%p where p is measured 
in GeV/c. Secondary tracks from decays and conversions 
are suppressed by matching tracks to hits in a second pad 
chamber (PC3) at distance of ~ 5.0 m. Track projections 
to the EMCal plane are used to veto photon candidates 
resulting from charged hadrons that shower in the EM- 
Cal. 



5 



III. METHOD 

A. Two-Particle Correlations 

Two-particle correlations are constructed by measur- 
ing the yield of particle pairs as a function of the mea- 
sured azimuthal angle between photon or parent me- 
son triggers and charged hadron partners. The corre- 
lation function, C(A0) = N^^^( (Ac^) / N^'^^^iAcf)) , cor- 
rects for the limited acceptance of 7-/1 or meson-hadron 
pairs by dividing the distribution in real events Nf^^^ 
by the mixed event distribution N^^" , . The corre- 

^ mixed 

lation function is decomposed utilizing a two-source 
model of pair yields coming from two-particle jet cor- 
relations superimposed on a combinatorial background 
yield from an underlying event. The underlying event 
in Au-|-Au is known to have an azimuthal asymmetry 
of harmonic shape quantified in the elliptic flow pa- 
rameter V2 [26l . [27j . This flow represents a harmonic 
modulation of the A(j) distribution of this underlying 
event, such that the flow-subtracted jet correlation sig- 
nal is encoded in the jet pair ratio function, JPR(A(j)) = 
C{A(j)) - £_{! + 2{vl){v^) COS2A0), using the notation of 
where (W2) is the average single-particle V2- 

Two methods of determining the background level ^, 
known as Zero- Yield at Minimum (ZYAM) and Abso- 
lute Normalization (ABS) respectively were applied to 
the Au-|-Au data. Both methods are described in detail 
in previous PHENIX pubhcations see also [i.[28l.[29| 
(ABS) and H (ZYAM). ZYAM assigns the level of zero 
jet yield and hence ^ to the minimum point of the corre- 
lation function C{A(f)). The ABS method uses the mean 
multiplicity of trigger-associated pairs in mixed events 
and a correction for finite centrality resolution to deter- 
mine ^. Where ZYAM statistical precision is reasonable, 
the direct 7-h extraction of the two methods agree to 
within much better than the total uncertainties, typi- 
cally within < 20%. The ABS method is chosen for the 
Au-|-Au results presented, as this method resulted in a 
more precise extraction of direct photon-jet pair yields 
at high trigger where lack of statistics near = 7r/2 
severely impairs the ZYAM determination. In the com- 
paratively low multiplicity p^-p collisions, the underly- 
ing event originates from different physical mechanisms 
than in Au-|-Au and is known not to be well described 
by event-mixing. Instead the correlation functions are 
normalized by fitting to a double Gaussian + constant 
function, corresponding to the ZYAM method [4] . 

The results presented here are corrected for the asso- 
ciated charged hadron efficiency eu such that the quoted 
yields correspond to a detector with full azimuthal ac- 
ceptance and I77I < 0.35 coverage. No correction is 
applied for the A77 acceptance of pairs. Final results 
are presented in terms of the yield Y of jet pairs per 
trigger, Y = A JPR{A(j))/Ntiiggcv with the constant 
^ = /<td(A</')/(27re.). ^ 

The magnitudes of elliptic flow were determined by 



measuring the distributions of inclusive photons, neutral 
pions, and charged hadrons as a function of the angle rel- 
ative to the reaction plane jwhich was determined with 
the BBC's as described in 31]. The V2 values measured 
for this analysis are consistent with previous PHENIX 
analyses m,!!^, H^]. At high-px ( > 6 GeV/c) the mea- 
sured tt'' V2 values used in the determination of the decay 
photon V2 are fit to a constant function in order to reduce 
the effects of large statistical fluctuations. The pt inde- 
pendence of V2 of TT^ 's is motivated by recent preliminary 
data [11] and also by the observed pr independence of 
the Raa, since parton energy loss is expected to be the 
dominant mechanism for V2 generation at high-px f^j . It 
is also consistent with the findings of [32] which is direct 
measurement of 7r° V2 for the same dataset and is being 
published concurrently with this measurement. Since, as 
discussed in that publication, the high-px functional be- 
havior for this dataset cannot be well-constrained, the 
level of uncertainty we assign to the constant fit assump- 
tion increases with px • 

Table H] lists the V2 values for the inclusive and tt" de- 
cay photons for all px ranges used, either the measure- 
ments, or for the highest px decay V2 values from the 
constant fit value. For the fit values the fit errors are 
listed as statistical error, despite the inherent systematic 
correlation of the fit value across the px bins. The de- 
cay photon V2 is derived from the measured tt'^ V2 by the 
same pl^ — > p^'^'^^'^ mapping procedure applied to the 
yields, described below. It is assumed that the V2 for 
other mesons which contribute decay photons {e.g. rf) 
are the same as that of the 7r° at high-px- This assump- 
tion is well motivated for the px range considered ( ^ 4.5 
GeV/c) under the expectation that the source of the high 
px azimuthal asymmetry V2 is jet quenching-induced sup- 
pression, already measured to be the same for a variety 
of mesons {e.g. r] itself [l^]) and by data measurements 
for other high px V2 which confirm the expectation [s^] 
for other hadrons. 



TABLE I: V2 values used in tlie jet function extraction for 
inclusive and decay photons in Au+Au collisions. 













Inelusive 


7 




Decay 7 


Centrality 


P 


'Y 
T 




Stat. 


Sys. 




Stat. 


Sys. 






5- 


-7 


0.053 


±0.009 


±0.011 


0.084 


±0.009 


±0.004 






7- 


-9 


0.047 


±0.022 


±0.015 


0.069 


±0.018 


±0.003 


0- 


20% 


9- 


12 


0.024 


±0.042 


±0.017 


0.069 


±0.020 


±0.003 






12- 


-15 


0.064 


±0.096 


±0.094 


0.069 


±0.023 


±0.003 






5- 


-7 


0.096 


±0.010 


±0.005 


0.155 


±0.011 


±0.036 






7- 


-9 


0.079 


±0.027 


±0.011 


0.105 


±0.019 


±0.025 


20- 


-40% 


9- 


12 


0.025 


±0.050 


±0.049 


0.105 


±0.020 


±0.025 






12 


-15 


0.287 


±0.128 


±0.104 


0.105 


±0.023 


±0.024 






5- 


-7 


0.143 


±0.023 


±0.035 


0.136 


±0.022 


±0.010 






7- 


-9 


0.146 


±0.064 


±0.026 


0.126 


±0.039 


±0.008 


40 


-60% 


9- 


12 


0.162 


±0.126 


±0.252 


0.126 


±0.042 


±0.008 






12 


-15 


-0.603 


±0.308 


±0.191 


0.126 


±0.046 


±0.008 
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B. Direct 7-Hadron Correlation Subtraction 

A direct photon is defined here to be any photon not 
from a decay process. Direct photons cannot be identified 
in Au+Au with reasonable purity on an event-by-event 
basis due to the large background of meson decay in the 
range of the analysis and the inability to use isolation 
cuts in the high multiplicity Au+Au environment. Thus 
both direct 7 and 7-/1 pairs must be determined from 
the already mentioned statistical subtraction procedure, 
which is therefore consistently used in this report for both 
the p+p and Au+Au. 

Single direct photons ha ve p reviously been measured 
in PHENIX, for Au+Au and p+p [H. In these 

analyses, the estimated yield of decay photons Nj^^^y is 
subtracted from a measured sample of inclusive photons 
-^inciusivo resulting in the direct photon yield. These mea- 
surements serve as an input to the current analysis, as 
they fix the fraction of the photon triggers which are ex- 
pected to be direct. This fraction is quantified by the 



fraction R~, 



inclusive 



/'^dccay ' 



The values used in 



this analysis are extracted from previous PHENIX mea- 
surements, [11,133] by interpolating to obtain the pt bin- 
ning used in this analysis. These interpolated values to- 
gether with the error estimations are tabulated in Table 

m 

TABLE II; Extracted R-^ values used as input to direct 7- 
h per-trigger yield subtraction (Equation [5}. These values 
are interpolated from previous PHENIX measurements as de- 
scribed in the text. 



Centrality 


Pt 




Stat. 


Sys. 




5-7 


1.77 


±0.09 


±0.06 




7-9 


2.45 


±0.09 


±0.18 


0-20% 


9-12 


2.99 


±0.11 


±0.41 




12-15 


3.66 


±0.24 


±0.68 




5-7 


1.46 


±0.10 


±0.04 




7-9 


1.85 


±0.10 


±0.12 


20-40% 


9-12 


2.30 


±0.12 


±0.28 




12-15 


2.35 


±0.20 


±0.44 




5-7 


1.30 


±0.09 


±0.05 




7-9 


1.52 


±0.07 


±0.13 


40-60% 


9-12 


1.85 


±0.10 


±0.30 




12-15 


1.94 


±0.24 


±0.36 




5-7 


1.18 


±0.01 


±0.06 




7-9 


1.33 


±0.01 


±0.05 


P+P 


9-12 


1.53 


±0.03 


±0.05 




12-15 


1.79 


±0.09 


±0.07 



The per-trigger yield of inclusive 7-/1 pairs ^inclusive is 
simply the weighted average of the contributions from 
decay and direct photon triggers. 



-^direct -^direct 



inclusive 



7 

decay-" decay 



Ide 



inclusive 



(1) 



Having already determined i?^, Idirect may then be ob- 
tained by simple manipulation of the above terms result- 
ing in statistical subtraction involving only per-trigger 
yields as follows. The decay photon per-trigger yield is 
subtracted from that of inclusive photons according to: 



R^Ys 



direct " 



'7 ^ inclusive 



-Y, 



decay 



R 



■7 



1 



(2) 



The direct 7 or direct 7-h pair yields do not, by defini- 
tion, exclude photons from jet fragmentation or medium 
induced photon production. 



C. Extraction of Decay Photon Correlations 

The decay photon associated yields are estimated from 
the measured Tr'^-h and 77-h correlations through a cal- 
culation which determines the decay correlations statis- 
tically from a Monte Carlo (MC) based, pair-by-pair 
weighting procedure. In this procedure the decay 7-h 
pair yield A^decay(PT) constructed by a weighted inte- 
gral over all Tr^-h and 77-h pairs. In what follows, we will 
first describe the procedure schematically, describing the 
ingredients and how they are obtained. We then give 
a more exact description and associated formula repre- 
senting exactly how the weighting was performed in the 
measurement. Schematically the procedure may be ex- 
pressed as a convolution of several factors according to 
the following relation, wherein for simplicity we only con- 
sider photons from tt*^ decay, although the procedure is 
also applied to rj decay photons. 



^t^{Pt) 



(3) 

where e^r and are the tt" and single decay photon ef- 
ficiencies, respectively, and V is the decay probability 
density, each of which is addressed in turn below. 

First, since the starting point is the uncorrected raw 
meson-h pair yield N'^^''', a correction for the parent me- 
son reconstruction efficiency, e^(p5')j is applied to the 
raw 7r° 's as a function of in order to account for the 
tt'^ daughter photons in the inclusive sample whose sis- 
ters lie outside the PHENIX acceptance or are otherwise 
undetected. Both efficiencies e-y, and e^^ in Equation [3] 
are also evaluated as a function of the position in the 
calorimeter along the beam direction, however this de- 
pendence mostly cancels in the ratio e^/e7r and therefore 
is suppressed for clarity. eT^ij)^) is determined by di- 
viding the raw number of tt'^'s N'^{p^) obtained in the 
same data sample by the PHENIX published tt" invariant 
yields @,[13,[3i| assuming no pseudorapidity dependence 
over the narrow PHENIX acceptance. The top panel 
in Fig. [1] illustrates, for the example of central Au+Au 
events, the 7r° efficiency correction factor l/eTr{p^). The 
correction rises at small px due to a px-dependent pair 
energy asymmetry cut designed to reduce combinatorial 
27 pairs reconstructed as real tt^'s. This cut, along with 
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FIG. 1: The weight factors used to obtain decay correlations 
from parent meson correlations. Top: vr'^ reconstruction effi- 
ciency correction, f/e,r. Middle: Decay probability function, 
Wab, for 5-7 GeV/c decay photons from vr" derived analyti- 
cally (black line), using the detector acceptance and resolution 
smearing (red line) and including the single decay photon effi- 
ciency, from a GEANT simulation (blue points). Bottom: 
^mcrgo Q^^ained by taking ratio of the blue points to red curve 
in the previous panel. 



the effects of any remaining background, is described be- 
low. At large px l/eir(PT) rises again due to losses from 
cluster merging. 

Second, the effect of decay kinematics is evaluated by 
determining the probabiUty density, V{p^,p'!^), for the 
decay of a px-independent distribution of tt^'s. 'P{p^,p^) 
represents the relative probabiUty of a tt" of p^—p^, to 
decay into a photon of p^. For a perfect detector, this 
function is calculable analytically. A simple fast MC 
generator implements the PHENfX acceptance and uses 
Gaussian smearing functions to simulate detector resolu- 
tion according to the known EMCal energy and position 
resolution. Occupancy effects give rise to an additional 
smearing of the tt'^ and ij invariant masses. This effect is 
included in the MC by tuning the resolution parameters 
to match the tt'^ peak widths observed in data. False re- 
construction of TT^'s and ?7's from combinatorial matches 
are either subtracted or assigned to the systematic un- 
certainties as discussed below. 

Finally, we wish to estimate the decay photon con- 
tribution to the measured raw inclusive photon sample 
which differs from the true decay photon distribution 
by the single decay photon efficiency, e^^p^). At in- 



termediate pt, e-y{p^) depends only on the photon mo- 
mentum and is included already implicitly by the fast 
MC simulation described above to produce V{p^,p^). 
Thus, it is useful to think of them as a single factor 

W^(Pt'Pt) = 'J-'{Ptt PT)^-y(PTT Pt) At high-PT, on the 
other hand, an efficiency loss is incurred by photons from 
tt'^'s whose showers merge into a single cluster in the 
calorimeter and are rejected by the shower-shape cut. As 
a consequence, the fraction of photons that are direct is 
artificially enhanced in the sample of reconstructed pho- 
ton clusters. The single decay photon efficiency depends 
on both the parent and daughter and is evaluated in a 
GEANT simulation, fn principle the convolution of both 
Vip^^p^) and €j{p^,p1^), W{p'^,p'!^), could be extracted 
as one function from the GEANT simulation, but ob- 
taining large enough MC statistics necessary to properly 
parameterize the above mentioned EMCal z position de- 
pendence of the e^^^ corrections is only feasible with the 
fast MC. Thus only the efhciency loss by cluster merging 
for photons e!^'"'s° is taken from the GEANT. The bot- 
tom panel of Fig. [T] shows e™°''8°(p^) evaluated from the 
GEANT simulation . 

Since we wish to construct per-trigger yields, the same 
procedure described in Equation [3] can be applied to find 
the estimated single decay photon trigger yield from the 
measured single tt^s, i.e. replacing N2~J^y with iV2ecay 
and N'^~^ with iV^. The exact application of schematic 
Equation [3] then takes the form of a sum over all 7r°-h 
pairs and single 7r°s found in the data. Each tt" or 7r°-h 
pair is given a weight which depends on tt*^ px- Opera- 
tionally we now split this weight into two parts: e7r(p^) 
discussed above and a factor Wabip^:)- The factor Wab is 
simply the end result of the fast MC-GEANT combined 
calculation, the convolution of V and e-^, including e™'^'^°, 
averaged over a chosen decay photon bin of the range 
a < pt < b. Thus in terms of the product W{p'^,p'!{) 
then Wafc(px) is given by 



Wab{PT) 



(4) 



Functions Wabip^^) are defined for the four photon px 
bins used in the analysis, [a,b] = [5,7], [7,9], [9,12] and 
[12,15] GeV/c. An example of Wabip^) for the 5-7 GeV/c 
bin is shown in Fig. [TJ Procedurally, we construct Wab as 
product of the fast MC curve shown in the middle panel 
and the linear fit discussed above to the bottom panel, 
gmcrgc j-pT -J ^ Although a decay of Px < the lower limit 
of the decay px bin, is kinematically disallowed, Wab is 
non-zero below this boundary when resolution effects are 
considered. Forp^ > 6, Wab decreases as ~ 1/Px ; slowly 
enough that tt^'s at values of px beyond the statistical 
reach of the data set contribute to the relevant decay pho- 
ton Px selections at a non-negligible rate. The tt^ sample 
is truncated at px — 17 GeV/c and extrapolated using 
power-law fits to the single and conditional tt" spectra 
to estimate a correction. In the latter case, each asso- 
ciated hadron px range is fit independently. The trun- 
cation avoids the high-px region where cluster merging 
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effects are dominant and the correction factor be- 
comes large. Although the truncation corrections for the 
number of decay photons and decay 7-/1 pairs are non- 
negligible, they mostly cancel in the per-trigger yield and 
are therefore typically < 1%, reaching a maximum value 
of 7% for only the 12 < pi^. < 15 » 3 < pJj, < 5 GeV/c 
bin. 

With the weight functions Wab the entire set of tt^- 
hadron pairs and single i:^ candidates (within a given 
range of A(/), < < (f)2, defining each Acj) bin) are 
then summed over, once for each decay photon px bin, 
and the per-trigger yield is constructed for each of these 
decay bins as 



i<pi<A4><4>2 



</-l<A0^_ft<02 



(5) 



In this form it is clear that the normalization of the func- 
tions eT(p^) and Wab{p^) cancel out completely in the 
per-trigger yield, and therefore only their shapes versus 
pl^ are important. Hence in Fig. [1] the curves are shown 
with arbitrary units. Also, as Equation [5] implies, the 
angular deviation between the direction of a decay pho- 
ton and its parent meson is ignored. The opening 
angle of a decay photon and hadron pair is taken to be 
the same as the A(l)^^h of the parent Tr^-h pairs. This 
approximation is tested in the fast MC and found to be 
extremely accurate since the distribution of angular devi- 
ation between a leading decay photon in a 27 decay and 
the parent mesons at these 7r° momenta have an RMS 
around of <C 0.01 radians, and the smallest A(p bins 
considered in the analysis are typically ~ 0.1 radians or 
larger. 



D. tt" and 77 Reconstruction 

In p+p collisions idccay is estimated using both recon- 
structed tt" and 77 mesons in invariant mass windows of 
120-160 and 530-580 MeV/c^ respectively. The total 
decay per-trigger yield is calculated from 



decay 



decay 



"/i/ttO-' decay 



(6) 



where <52/7rO ratio of the total number of decay pho- 

tons to the number of decay photons from vr*'. Based on 
the measurements of 77 [2^ and uj [s^, which together 
with the tt" account for > 99% of decay photons, the 
value of <5^y^o is determined to be 1.24±0.05 in the high- 
Pt region covered by this analysis, independent of col- 
lision system and centrality. Note that the per-trigger 
yields for to and other heavier meson triggers (uj,rif,(j),...) 
are not measured and are taken to be equivalent to Y^^ 
in Equation [HI This assumption was studied in PYTHlA 
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5-7 ® 2-3 GeV/c 
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FIG. 2: (color online) Examples of parent and daughter per- 
trigger yields for the tt" and rj in p+p collisions for px se- 
lection 5 < 7 and 2 <px < 3 GeV/c. These correlation 
measurements are used to determine the total decay photon 
per-trigger yield as described in the text. 



and found to influence Fdocay at the level of < 2%. In 
Au-|-Au collisions correlations using 77 triggers are not 
directly measured, but rather estimated from the p+p 
measurement as discussed below. 

Figure [2] shows the various components of the de- 
cay photon measurement in p+p. In p+p collisions the 
rate of combinatorial background photon pairs is re- 
duced by only considering photons of px > 1 GeV/c 
resulting in background levels of < 10% for which no 
correction was applied. The effect of such remaining 
pairs on i^J^^ay evaluated to be negligible (< 2%) 
compared to the size of other uncertainties on the final 
^direct result using a detailed full PYTHIA test of the 
method which included 7r° reconstruction with combina- 
torial photon pairs. On the other hand, 77 reconstruc- 
tion has a much smaller signal-to-background of 1.4-1.6, 
depending on the px selection, even in the low multi- 
plicity p+p environment. In this case, the per-trigger 
yield of the combinatorial photon pairs is estimated from 
photon pairs with invariant mass in "sideband" ranges 
of 400-460 and 640-700 MeV/c^ beyond 3cr of the 77 
peak. The sideband contribution i^^ecay'^"'^ then sub- 
tracted using the signal-to-background ratio /bkg eval- 
uated from gaussian -I- polynomial background fits to 
the invariant mass distributions according to ^g^ay' ~ 

^^JcTay/lV/bkg + 1) - lli^:^^"V/bkg. The yield Fdtay""' 

is generated from the full meson to decay photon weight- 
ing function procedure (Equation [5]). The subtraction 
procedure was also tested in PYTHIA and the extracted 
and input per-trigger yields were found to agree to within 
10%. 

In Au-|-Au collisions the combinatorial rate for tt" re- 
construction is substantially larger. Correspondingly, a 
Px dependent cut on the pair energy asymmetry ai2 = 
|£;i -£;2|/(£;i -^£^2) [23II, visible in Fig. [U with the small- 
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est allowed asymmetry at the lowest 7r° pT values, is used 
to reduce this background. With such cuts the signal-to- 
background in central events varies from 5:1 at its lowest, 
increasing to about 15:1 for the highest pT selection. The 
effect of the combinatorial background is studied through 
examination of a similar sideband subtraction analysis as 
in the p+p rj — h correlation extraction described, this 
time for Tr^-h, using invariant mass ranges just outside 
the tt" peak region. However no clear trend beyond non- 
negligible statistical limitations is observed, so no correc- 
tion for the background is applied. Instead the maximum 
size of the effect (typically ~ 7%) is included as source 
of systematic uncertainty on the decay yields and prop- 
agated to the final direct photon per-trigger yields. 

In central Au+Au collisions the rj meson cannot be 
reconstructed with sufficient purity to measure its corre- 
lations. Instead, a scaling argument is employed. Moti- 
vated by the similar high-px suppression pattern shown 
by rj and tt" in Au-|-Au [24] and corresponding near 
equality of the p+p and Au-|-Au ij/tt'^ ratios, the ratio 
lCy(^)/lCy(^o) is measured in p+p and applied as a cor- 
rection to the Au-|-Au 1^(7^"). This is justified by the 
assumption that the jet fragmentation is primarily oc- 
curring outside the medium. We do not attribute any 
additional uncertainty to this scaling beyond the 10% 
sideband systematic and statistical uncertainties of the 
r] measurement in p+p. However, to give an idea of the 
possible impact of this assumption, the total systematic 
error on Idccay from all other sources would correspond 
to a variation of the Au-|-Au by ~ 50%. Given the 
similarity of the high-pT suppression demonstrated by all 
light quark bound states measured thus far, this would 
correspond to a rather large change. 



IV. SYSTEMATIC UNCERTAINTIES 

There are four main classes of systematic uncertainty 
in the Au-|-Au data: elliptic fiow, normalization of the 
underlying event (ABS), R-y, and the decay per-trigger 
yield estimate, the latter two of which are present in the 
p+p data as well. Table IIIII lists the fractional contri- 
bution of each of these sources to the total systematic 
uncertainty on the direct photon per-trigger yields in the 
20% most central Au-|-Au and p+p data. In the cen- 
tral Au-|-Au data the uncertainty at low p!^ is dominated 
by the V2 and correlation function normalization (ABS 
method) estimation due to large multiplicity of hadrons. 
At higher pi^, but low trigger px, P^t the decay error dom- 
inates due to the two-photon combinatorial background 
for 7r° reconstruction. Finally, at large pt^ and ^he 
backgrounds responsible for both of these sources of un- 
certainty decrease and the uncertainty on i?^, which is 
relatively constant, dominates. In p+p collisions the de- 
cay photon background forms a much larger fraction of 
the total photon sample. In this case, the decay uncer- 
tainty arises from the MC decay photon mapping proce- 
dure, the ry sideband subtraction and the rj/n'^ ratio in 



TABLE III: Fractional contribution to the total systematic 
uncertainty for each of the main sources of uncertainty in 
p+p and 0-20% Au+Au collisions. 





Pt 


Au-(-Au, Ccntrality 0-20 % 


P+P 




(GeV) 


(GoV) 




Decay 


v2 


Norm. 


R-, 


Decay 




1-2 


0.03 


0.14 


0.50 


0.33 


0.14 


0.86 


5-7 


2-3 


0.02 


0.32 


0.46 


0.20 


0.21 


0.79 




3-5 


0.02 


0.71 


0.18 


0.10 


0.05 


0.95 




1-2 


0.09 


0.17 


0.45 


0.29 


0.22 


0.78 


7-9 


2-3 


0.10 


0.35 


0.38 


0.17 


0.25 


0.75 




3-5 


0.09 


0.61 


0.18 


0.13 


0.21 


0.79 




1-2 


0.06 


0.09 


0.53 


0.33 


0.19 


0.81 


9-12 


2-3 


0.26 


0.25 


0.33 


0.16 


0.30 


0.70 




3-5 


0.46 


0.30 


0.13 


0.10 


0.35 


0.65 




1-2 


0.08 


0.01 


0.63 


0.29 


0.21 


0.79 


12-15 


2-3 


0.21 


0.14 


0.48 


0.17 


0.02 


0.98 




3-5 


0.22 


0.14 


0.39 


0.25 


0.10 


0.90 



approximately equal parts. The yields associated with 
daughter photons are larger than for the meson parents 
because of feed-down from larger values of parent px, and 
hence, jet px- 

The correction for single hadron efficiency e;i(px) 
varies as a function of collision system and centrality. 
These corrections are obtained by finding the ratio of 
raw yields of hadrons obtained without the trigger con- 
dition in the same analysis (i.e.) with the same cuts as 
in the analysis, to the previous PHENIX published mea- 
surements of the corresponding charged hadron spectra. 
[40I [ilj . As in previous PHENIX two-particle correla- 
tion measurements, 0, |30|, this procedure has inherent 
uncertainties assigned as a px-independent 10% uncer- 
tainty, on each system and/or centrality. 

V. RESULTS 
A. Direct 7-h Per-Trigger Yields 

Figure [3] shows examples of direct photon per-trigger 
yields in p+p and central Au+Au collisions. Also shown 
are the per-trigger yields for inclusive and decay pho- 
ton triggers which are the ingredients in the statistical 
subtraction method as expressed in Equation [21 A clear 
away-side correlation is observed (Acji ~ tt) for direct 
photons triggers in p+p. In Au+Au collisions the away- 
side correlation is suppressed for both decay and direct 
photon triggers. The near-side direct photon associated 
yields are small relative to that of decay photons, an ex- 
pected signature of prompt photon production [16]. 

The away-side yields, integrated over \A<j) — tt] < 7r/5 
radians, are shown in Fig. [3] and Table IIVI for p+p and 
Au+Au collisions. This range roughly corresponds to 
the "head region" as defined in Q and is chosen primar- 
ily to minimize the influence of medium response which 
is thought to dominate the "shoulder" region further off- 
set from A(j> = IT. Additionally, the acceptance and the 



10 



0.2 

^ 0.1 

0^ 



T3 



-e- 
< 

T3 



0.2 
0.1 




p+p 



Inclusive y-h 
□ Decay y-h 
• Direct y-li 



Au+Au, 0-20% 



□ 

u> 



0.5 1 1.5 2 
A(^ [rad] 



2.5 



FIG. 3: (color online) Examples of per-trigger yields used in 
the direct photon correlation analysis for the 5 < 7 and 
3 <Pt < 5 GeV/c bin. Top (bottom) panel: Inclusive, decay 
and direct photon per-trigger yields in p+p (0-20% central 
Au+Au) collisions. 
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FIG. 4: (color online) Direct 7-/1 per-trigger yields for the 
range \A<I) — n\ < 7r/5 radians vs. associated hadron pt- 
Four different direct 7 px ranges (indicated on the figure) are 
shown in the most central 20% of Au-|-Au events and p+p 
events. The upper limits are for 90% confidence levels. A px- 
independent uncertainty of 10% due to the charged hadron 
efficiency correction is suppressed from the plot. 



TABLE IV: Direct -y-h per-trigger yields in 20% most central 
Au+Au and in p+p collisions. An additional px-independent 
uncertainty of 10% due to the charged hadron efficiency cor- 
rections is not shown. 



ph 



Pt 

(GeV) (GeV) 



Yield 
Au+Au, 



Stat Sys 
Centrality 0-20% 



Total 



5-7 



1- 2 0.23 

2- 3 0.41 

3- 5 0.62 



6.26e-02 4, 
2.68e-02 1, 
4.82e-03 2, 



72e-02 
29e-02 
13e-03 



4.62e-02 
5.68e-03 
1.96e-03 



6.60e-02 
1.41e-02 
2.90e-03 



7-9 



1- 2 0.17 

2- 3 0.3 

3- 5 0.46 



3.71e-02 8, 
3.45e-02 2, 
9.63e-03 4, 



48e-02 
39e-02 
18e-03 



5.59e-02 
8.46e-03 
1.96e-03 



1.02e-01 
2.53e-02 
4.62e-03 



9-12 



1- 2 0.13 

2- 3 0.23 

3- 5 0.36 



1.28e-01 1, 
3.94e-02 3, 
-2.16e-03 6, 



34e-01 
81e-02 
29e-03 



6.84e-02 
l.Ole-02 
2.06e-03 



1.51e-01 
3.94e-02 
6.62e-03 



12-15 



1- 2 0.1 

2- 3 0.18 

3- 5 0.28 



5.31e-01 2, 
-6.13e-03 6, 
3.25e-02 1, 



53e-01 
99e-02 
60e-02 



1.49e-01 
1.80e-02 
2.52e-03 



2.94e-01 
7.22e-02 
1.62e-02 



p+p 



5-7 



1- 2 0.24 

2- 3 0.43 

3- 5 0.66 



1.44e-01 9, 
4.22e-02 5, 
1.55e-02 2, 



93e-03 
47e-03 
07e-03 



3.42e-02 
1.20e-02 
3.26e-03 



3.56e-02 
1.32e-02 
3.86e-03 



7-9 



1- 2 0.18 

2- 3 0.31 

3- 5 0.48 



1.73e-01 1, 
6.24e-02 1, 
2.26e-02 4, 



84e-02 
lle-02 
53e-03 



2.88e-02 
1.15e-02 
3.75e-03 



3.42e-02 
1.60e-02 
5.88e-03 



9-12 



1- 2 0.14 

2- 3 0.24 

3- 5 0.38 



2.59e-01 2. 
7.01e-02 1. 
1.94e-02 7. 



99e-02 
73e-02 
21e-03 



2.50e-02 
l.OOe-02 
3.77e-03 



3.90e-02 
2.00e-02 
8.14e-03 



1-2 0.11 
12-15 2-3 0.19 
3-5 0.3 



1.20e-01 5, 
1.04e-01 3, 
4.26e-02 1, 



13e-02 
lle-02 
62e-02 



7.22e-02 
2.02e-02 
1.13e-02 



8.86e-02 
3.71e-02 
1.97e-02 



the head region. Due to statistical and systematic fluctu- 
ations, the subtraction of the decay-photon hadron pairs 
from the inclusive 7-/1 sample can result in a negative 
yield. In this case 90% confidence-level upper limits are 
given. In the case that a positive yield is obtained, but 
the uncertainty is consistent with 0, the lower bound of 
the error bar is also replaced with an arrow. As noted 
in the figure caption, a 10% px-independent uncertainty 
due to the charged hadron efficiency corrections is not 
shown. 



signal itself are largest in this range so statistical preci- 
sion is maximized. It should be noted that the width of 
the jet correlation is larger than this interval. We do not 
make a correction for this effect, since we are primarily 
concerned with the comparison of the yields from p+p 
and Au+Au collisions. It should be noted, however, that 
in addition to parton energy loss, any broadening of az- 
imuthal correlations, whether by hot or cold nuclear mat- 
ter effects, will contribute to a suppression in the yield in 



B. Suppression Factor Iaa 

Departure from the vacuum QCD processes is quan- 
tified by Iaa, the ratio of Au+Au to p+p per-trigger 
yields: 

/aaW,Pt)- y,+p(p7^p.) ■ 
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FIG. 5: (color online) Ratio Iaa of the Au+Au to p+p yields 
shown in Fig. |4] An additional pT-independent uncertainty 
of 14% due to the charged hadron efficiency corrections is not 
shown. 
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FIG. 6: (color online) /aa(Pt ) integrated over the range 5 < 
p^ < 15 GeV/c for associated hadrons of 3 < < 5 GeV/c 
vs. centrality compared to single tt" high-px Raa (integrated 
overpT > 5 GeV/c) An additional pT-independent uncer- 
tainty of 14% due to the charged hadron efficiency corrections 
is not shown. 



Figure [5] shows the Iaa values for all direct photon and 
associated hadron bins for the most central 0-20% of 
collisions. The data points for which the subtraction 
resulted in a negative yield value (the 90% confidence 
level upper limits) are included with standard 1-a un- 
certainties. For the range 5-12 GeV/c, a significant 
suppression is observed in the 3 < pt^ < 5 GeV/c bin 
in which the highest precision is obtained. At lower p^, 
where the background subtraction is largest, the data do 
not have the statistical precision to determine the degree 
to which the yields are suppressed. Iaa for direct photon 
triggers is consistent to that of charged hadron triggers 
[3| as shown in the top left panel in which results with 
similar ranges of pT,t are compared. 

Figure [S] shows the Iaa for the pJ^ = 3-5 GeV/c bin, 
integrated for all trigger pT bins {p^ = 5-15 GeV/c) 
and for three centrality bins, 0-20%, 20-40%, and 40- 
60%. For the most central bin, the suppression of the 
away-side direct photon per-trigger yield is clearly ob- 
served, Iaa = 0.32 ± 0.12'^''^* ± 0.09*^^"*. Within large 
uncertainties we see that the 7-jet Iaa in this pt range, 
dominated by moderate to high values of z (= p''/p'°'), 
is consistent with the single particle Raa as a function of 
centrality, consistent with a scenario in which the geom- 
etry of suppression plays an important role as would be 
expected from a sample dominated by surface emission. 

Figure [6] also compares Iaa from a measurement of 
high-pT dihadron {h^ — h^) correlations Q to the 7- 
jet result for similar px^t selections. The two results are 
remarkably similar in the most central bin. This may in- 
dicate that surface emission is dominant for both samples 
in this z region. However it should be noted that the total 
uncertainties on either measurement are still quite large 
on a relative scale. As explained in the introduction, the 



two measurements should be subject to different geomet- 
rical effects. Disentangling such effects through precise 
comparisons of dihadron and 7-/1 suppression should be 
pursued with future measurements with improved statis- 
tics. 



C. Towards the Fragmentation Function 

Using the distribution of charged hadrons opposite di- 
rect 7 triggers, parton energy loss may be studied di- 
rectly as a departure from the (vacuum) fragmentation 
function. In distinction to Tr'^-h correlations, where the 
away-side distribution is only sensitive to the integral of 
the fragmentation function (the average multiplicity of 
the away-side jet) [i^] , the away-side distribution for di- 
rect 7-h correlations provides a measurement of the full 
fragmentation function of the jet from the away-side par- 
ton. To the extent that the transverse momentum of the 
away-side parton and the direct 7 are equal and opposite, 
as in leading order pQCD, the fragmentation function of 
the jet from the away-parton should be given to a good 
approximation by the xe distribution, 

-Pt ■ Pt ~Pt cos Ac/) 

where the transverse momentum of the trigger pt^, = p^ 
in the case of 7-h correlations. The reasons why the scal- 
ing variable xe is an approximation to, rather than exact 
measure of, the fragmentation variable of the away-side 
jet with momentum Za are: i) the away-side parton does 
not generally balance longitudinal momentum with the 
trigger 7, although it is restricted by the Ary acceptance 
of the detector; ii) the transverse momenta of the 7 and 
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FIG. 7: zt distributions dN/dzT from the direct photon as- 
sociated yields in p+p (left) and 0-20% Au+Au (right) colli- 
sions. 

away parton do not exactly balance. The transverse mo- 
mentum imbalance was discovered at the CERN-ISR us- 
ing xe distributions |43j| and originally attributed to an 
"intrinsic" transverse momentum /ct of each of the initial 
colliding partons [i^l, but now understood to be due to 
"resummation" of soft-gluon effects [i^ . 

The validity of the approximation xe ^ Za can be 
tested by observing identical xe distributions for dif- 
ferent values of trigger {xe scaling), in which case 
one would accept the xe distribution in 7-h correlations 
as the quark fragmentation function from the reaction 
g + (?—!■ g + 7 without need of correction. We approxi- 
mate Xe zt, the ratio of the mean associated to 
mean trigger for each hin.^ The (p^) for the four 
trigger bins are: 5.66, 7.75, 10.07, 13.07 GeV/c, close to 
the values obtained from a fit to the direct-7 invariant 
cross section of the form p^^^ j35l] . 

Figure [7| shows the zt distributions for p+p and 
Au-f Au collisions. The p+p data (Fig. [7^) exhibit rea- 
sonable Zt scaling so that the measured distribution 
should represent the away-side jet fragmentation func- 
tion. A fit of this data to a simple exponential {Ne~''^'^) 
gives an acceptable y^/dof = 12.8/10 with a value 
h = 6.9 ±0.8 which is consistent with the quark fragmen- 
tation function, parameterized [42| as a simple exponen- 
tial with 6 = 8.2 for 0.2 < z < 1.0, and inconsistent with 
the gluon fragmentation function value of b ~ 11.4. It 
should, however, be recalled that the data do not cover 
the full extent of the away peak, only |A0 — 7r| < 7r/5 
radians, and that possible variations of the widths of 
the peaks in both the p+p data and the Au-|-Au data 



^ The reader is advised to carefully distinguish this variable = 
{pt^)/{pij,) from our previous notation used in [4^ of zt = ptj,/p, 
which is the fraction of jet momentum p contained in the trigger 
particle. 



FIG. 8: (color online) Iaa{zt) for the 20 % most central 
Au+Au data compared to predictions from an energy loss 
calculation An additional pT-independent uncertainty 

of 14% due to the charged hadron efficiency corrections is not 
shown. 

with p^ and pi^ have not been taken into account in the 
present analysis. Additionally a more detailed analysis, 
differential in trigger p^ , is necessary to study trigger p^ 
dependent effects which can influence the fragmentation 
function fit values [43 |. 

In central Au+Au collisions, the fragmentation func- 
tion may be modified by the medium^ , so that Zt scaling 
should not hold except in two special cases: i) pure sur- 
face emission or punch-through where the away-side jets 
are not modified — the zt distribution will be suppressed, 
but will have the same shape as in p+p collisions; ii) 
constant fractional energy loss of the away jet — the zt 
scaling will be preserved in Au+Au collisions but with a 
steeper slope than in p+p collisions. The Au+Au data 
(Fig. [Tb) are consistent with zt scaling with the same 
shape as the p+p data, a value of 6 = 5.6 ± 2.2 and ex- 
cellent x^/dof = 10.1/10 for the simple exponential fit. 
The point at lowest zt = 0.11 for Au+Au is 1.6 stan- 
dard deviations above the fit, suggesting that improved 
statistics will permit the observation of any non-surface 
emission. 



D. Model Comparison 

Several authors have reported predictions for 7-jet in 
heavy ion collisions [i^, H^, [13, [EH . As a demonstration 
of the how such calculations can be compared to the data, 
the /aa values as a function of zt are compared to energy 
loss predictions 49] in Fig. [51 The calculation uses effec- 
tive fragmentation functions to parameterize the energy 



2 See Equation 1 in [i^l 
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loss in terms of a parameter Eq which is expected to be 
proportional to the initial gluon density [131 • The model 
calculates the energy-loss of the leading parton, and ne- 
glects the contribution the gluon radiation and medium 
response which may dominate at low values of z. The 
data is well reproduced by the model over the range of 
values of Sq provided, 1.48-1.88 GeV/fm. This corre- 
sponds roughly to the range of Sq allowed by comparison 
to the PHENIX Raa data of 1.9to;5 

It should be noted that the calculation rejects frag- 
mentation photons with an isolation cut. Such a proce- 
dure has not yet been demonstrated in central Au-|-Au 
data, although doing so would help to eliminate beyond- 
leading-order effects. 

VI. CONCLUSIONS 

We have presented the first direct 7-h measurements 
in Au+Au and p+p collisions at RHIC. A significant sup- 
pression of Iaa = 0.32 ± 0.12'^'^* ± 0.09"^"* for the away- 
side charged hadron yield in the range 3 < < 5 GeV/c 
is observed for direct photon triggers in Au-|-Au as com- 
pared to p+p. Furthermore, the level of suppression is 
found to be consistent with the single particle suppres- 
sion rate and the importance of energy-loss geometry, 
notably the expectation of surface emission in the kine- 
matic range sampled. A possible indication that energy- 
loss geometry may also be important in dijet suppression 
is that "f-h suppression Iaa is also observed to be quite 
similar to that of dihadron suppression in central events; 
however, the current precision of the data does not ex- 
clude substantial differences. In the p+p data zt scaling 
is observed, suggesting that the measured zt distribution 
(Fig. [7]) is a statistically acceptable representation of the 
fragmentation function of the quark jet recoiling away 
from the direct photon. Improvement of the statistical 
and systematic precision of the measurements should al- 
low further tests of vacuum fragmentation expectations 
in p+p collisions and insights into details of the medium 
modification of jet fragmentation in Au-|-Au. 
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